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We analyze data for ^^C(e, e'p) with < 2 {GeY/cf using the relativistic distorted- wave impulse 
approximation (RDWIA) based upon Dirac-Hartree wave functions. The Ip normalization extracted 
from data for > 0.6 (GeV/c)^ is approximately 0.87, independent of Q^, which is consistent 
with the predicted depletion of the IP3/2 orbital by short-range correlations. The total Ip and Is 
strength for Em < 80 MeV approaches 100% of IPSM, consistent with a continuum contribution 
for 30 < < 80 MeV of about 12% of IPSM. Similarly, a scale factor of 1.12 brings RDWIA 
calculations into good agreement with ^^C(e,e'p) data for transparency. We also analyzed low 
data from which a recent NRDWIA analysis suggested that spectroscopic factors might depend 
strongly upon the resolution of the probe. We find that the momentum distributions for their 
empirical Woods-Saxon wave functions fit to low data for parallel kinematics are too narrow to 
reproduce data for quasiperpendicular kinematics, especially for larger Q^, and are partly responsible 
for reducing fitted normalization factors. Although the RDWIA normalization factors for < 0.2 
(GeV/c)^ are also smaller than obtained for > 0.6 (GeV/c)^, the effect is smaller and we argue 
that it should be attributed to the effective single- nucleon current operator instead of to spectroscopic 
factors which are probe-independent properties of nuclear structure. However, remediation of the 
failure of RDWIA calculations to reproduce low data for parallel kinematics will require a more 
sophisticated modification of the current than a simple multiplicative factor. 



I. INTRODUCTION 



It is generally believed that single-nucleon electromagnetic knockout provides unambiguous measurements 
of the momentum distributions and spectroscopic factors for single-hole states near the Fermi surface. Re- 
views of analyses based upon the nonrelativistic distorted-wave impulse approximation (NRDWIA) can be 
found in Refs. which show that the momentum distributions are described well by mean-field cal- 

culations, such as Skyrme-Hartree-Fock, while spectroscopic factors for low- lying state are reduced relative 
to the independent particle shell model (IPSM) by an average factor of about 65% over a broad range of A. 
Part of the depletion of IPSM strength can be attributed to short-range correlations that shift approximately 
15% of the hole strength to very large missing energies, beyond 100 MeV. Recently, a direct measurement 
of the spectral function using ^^C(e,e'p) in parallel kinematics ;4] observed approximately 0.6 protons in a 
region with pm ^ 240 MeV/c and Em ^ 50 MeV attributable to single-nucleon knockout from correlated 
clusters. This finding is consistent with the 16% IPSM predicted by Frick et al. Q using the self-consistent 
Green's function (SCGF) method Q. Similar estimates of the depletion of hole states are available from 
the correlated basis function (CBF) theory and other methods also. Furthermore, SCGF calculations 
1^ 01 show that the momentum distribution for low-lying quasihole states is hardly affected by short-range 
correlations, remaining very similar to mean-field calculations, while the momentum distributions for the 
large- £',„ continuum are substantially broadened. 

If the occupancy of IPSM orbitals is predicted to be approximately 85% by theories that reproduce the 
correlated strength at large {Em,Pm), why is only 65% observed by (e, e'p) experiments? Coupling to collec- 
tive modes is expected to produce significant fragmentation of the valence quasi-hole strength spread over a 
range of perhaps 10 MeV, but the extended RPA calculations presently available do not yet reproduce those 
fragmentation patterns well @. Thus, it is possible that many small fragments are missed experimentally 
even if they lie within the experimental acceptance. Alternatively, problems in the reaction model may lead 
to systematic underestimation of spectroscopic factors. Most of the data contributing to the aforementioned 
estimate of 65% IPSM were limited to Tp < 100 MeV and were analyzed using NRDWIA calculations based 
upon empirical Woods-Saxon wave functions. More recent relativistic analyses (RDWIA) typically produce 
larger spectroscopic factors and describe interference response functions more accurately 0, 0, . 
The effects due to dynamical enhancement of Dirac spinors by the nuclear mean field are actually stronger 
at low than at higher energies [T^ . 

While spectroscopic factors are of obvious importance to theories of nuclear structure, they also affect 
analyses of nuclear transparency that look for the onset of color transparency 0, . Those experiments 
typically compare the (e,e'p) yield integrated over wide but finite ranges of missing momentum, pm, and 
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missing energy, E„i, with calculations based upon a model spectral function. This spectral function is usually 
based upon the IPSM with a correction factor for the depletion by correlations of the single-hole strength 
within the experimental acceptance. Therefore, the accuracy of the depletion factor affects the accuracy of 
transparency measurements. 

Recently Lapikas et al. Tt'I analyzed the data for i2c(e,e'p) with < 0.4 {GeV/cf using NRDWIA 
calculations with Woods-Saxon bound-state wave functions fit to the NIKHEF data for Tp — 70 MeV in 
parallel kinematics |0| and concluded that the normalization factors relative to IPSM were Nip — 0.56 ±0.02 
and Nis = 0.59 ± 0.04 for Em < 80 MeV. They then argued that the transparency results based upon an 
estimated depletion factor of / = 0.9 for the single-hole spectral function with Em < 80 MeV should be 
increased by a factor of approximately 1.25 using appropriately weighted attenuation factors for the Ip and 
Is shells. Finally, they argued that the normalization factors for ^^C vary strongly with from the low 
values fitted to the NIKHEF data to values approaching the IPSM limit for > 2 (GeV/c)^. Frankfurt 
et al. ^3 speculated that quenching of spectroscopic factors for low might be caused by probing a 
dressed quasiparticle using long wave lengths and that at higher resolution the effective current operator 
for the (e, e'p) reaction approaches more closely that for a free nucleon. However, although the range 
was smaller, our recent relativistic analysis showed no significant slope in the normalization factors fit to 
1^0(6, e'p) data for < 0.8 (GeV/c)^ Furthermore, wc find that RDWIA calculations based upon the 
ordinary single-nucleon current with dynamical enhancement of lower components of Dirac spinors describe 
the left-right asymmetry for quasiperpendicular kinematics very well [T^ : it is not obvious that one should 
describe quasiparticles for low using Dirac spinors in a similar mean field. 

In this paper we use relativistic distorted wave (RDWIA) calculations to analyze the data for ^^C(e, e'p) 
with < 2 (GeV/c)^. In addition to the data considered by Lapikas et al. [ij, we include recent data 
from Dutta et al. ^3 fo'^ — 0-6, 1.2, and 1.8 (GeV/c)^. Our RDWIA fits produce systematically larger 
normalization factors using Dirac-Hartree wave functions than using the Woods-Saxon wave functions that 
were fit to the NIKHEF data. The latter do not describe the momentum distributions for larger well; nor 
do they describe the data for low with quasiperpendicular kinematics as well as do Dirac-Hartree wave 
functions. Therefore, the fits to data for parallel kinematics with small ejectile energies appear to produce 
artifically narrow momentum distributions that result in anomalously low spectroscopic factors. 

The model and fitting procedures are described in Sec. ^ Sec. IIIII presents an analysis based upon 
the effective momentum approximation (EMA) that permits direct comparison between nonrelativistic and 
relativistic bound-state wave functions. Sec. IIVI presents a more rigorous analysis, without using the EMA, 
of both spectroscopic factors and transparency. Our conclusions are summarized in Sec. 

II. MODEL AND DATA 

Detailed descriptions of RDWIA for (e, e'p) reactions are been given in many recent papers. We will use the 
methods and terminology of Refs. jT^,^^!. All calculations treat electron distortion in the qcs approximation 
and use the CC2 current operator with MMD form factors in Coulom b g auge. Optical potentials for 
Dirac phenomenology were taken from the global analysis by Cooper et al. [22| ■ We consider Dirac-Hartree 
wave functions from the original Horowitz and Serot (HS) analysis '23*1 and NLSH wave functions from 
Sharma et al. (2J|]. In addition we consider Woods-Saxon wave functions fit to ^^C(e, e'p) data. 

The data we consider are summarized in Table |l| Reduced cross sections for the Tokyo, Saclay, and 
SLAG data were provided by Lapikas T^, '25j and include small adjustments to a common convention for 
reduced cross section. In addition, radiative corrections were applied to the data from Ref. 26] . Only for 
the NIKHEF experiment |0| was it possible to resolve the lowest three Ip fragments. Although there are 
small differences between the missing momentum distributions for these three fragments, the ground state 
represents about 81% and the next lp3/2 fragment an additional 9% of the total Ip strength for Em < 25 
MeV. Therefore, we represent the Ip strength using either a single Dirac-Hartree wave function for IP3/2 or 
the NIKHEF fit to the ground state. That experiment was also able to resolve several weak positive-parity 
states with Em < 30 MeV, but their strength is small enough to neglect in the present analysis (23 ■ For 
most other experiments the lowest missing energy bin contains a small contribution from the low-energy 
tail of the Is shell. Furthermore, for all experiments the bin intended to emphasize the Is shell inevitably 
contains contributions from a broad continuum that may include additional Ip strength. Therefore, each 
Em bin will be fit as an incoherent mixture of Ip and Is contributions according to 



Crcd = A^lpO-rcd(lP3/2) + A^lsCTrcd (lSi/2) 



(1) 
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where the overlap function for each shell is assumed to be independent of missing energy and is normalized 
to full occupancy. In addition, for the JLab data sets we fit a normalization factor for the independent 
particle shell model (IPSM) to the inclusive data in the range 10 < Em < 80 MeV according to 

CTrcd = ^IPSM (o-rcd(lj3) + CTred(ls)) (2) 

Calculations will be shown for the entire Pm range measured for each data set, but only data in the range 
\Pm I < 200 MeVc will used to fit the normalization factors because larger pm may be susceptible to increasing 
corrections for effects not contained in the RDWIA model based upon single-nucleon current operators. 

We also require a model of the Is energy distribution to obtain the spectroscopic factor from the nor- 
malization A'^is(i?min, £'max) fit to limited and variable intervals of Em- One often employs a Lorentzian 
distribution 

L,s{Em) = ^ ^^^^ _ ^J, +"(r(i?,„)/2)2 

whose energy-dependent width is given by the Brown-Rho parametrization 

[Era - EfY + b 

where Ep is the Fermi energy and a,b are constants. For the present purposes it suffices to use the parameters 
a = 24 MeV and b = 500 MeV^ originally proposed by Brown and Rho "^^l with Ep = 15.96 MeV and 
Els = 39 MeV for ^^C. Lapikas et al. obtained slightly different parameters for b and Eis by fitting the 
Tokyo data, but the effect upon Sis is insignificant. Thus, one would expect approximately 11% of the 
Is strength to be broadly distributed above 80 MeV missing energy, which is consistent with the depletion 
factor attributed to short-range correlations. Therefore, in the absence of background, the Is spectroscopic 
factor deduced from a fitted normalization factor would be 

\ Jo ^ls{Em^dEm 

Sis = 2A^ls(£^min, £^max) E^^^^ ("^^ 

J E """^ Lis{Em)dEm, 



TABLE I; Summary of data for ^^C(e,e'p). 
Quasiperpendicular kinematics with constant 



Eo is the beam energy and Tp is the central proton kinetic energy. 
{uj, q) were used unless a range is given. 



data set 


Eo 


AS™(lp) 


A£;„(ls) 




Q2 


notes 




MeV 


MeV 


MeV 


MeV 


(GeV/c)2 




Tokyo [29] 


700 


6-30 


21-66 


159 


0.29 






700 




21-66 


136 






Saclay gO] 


497 


15-22 


30-50 


87 


0.16 




Saclay [26] 


500 


15-22 


NA 


99 


0.18 


constant {to, q) 


Saclay ]26] 


500 


15-22 


NA 


99 


0.09-0.32 


parallel; Is not available 


NIKHEF [18j 


285-481 


15-22 


30-39 


70 


0.02-0.26 


parallel 


SLAC [^Ij 


2015 


15-25 


30-80 


600 


1.11 


Is not radiatively unfolded 


JLab [15j 


2445 


15-25 


30-50 


350 


0.64 


inclusive bin with 10 < Em < 80 MeV available 


2445 






700 


1.28 






3245 






970 


1.84 





III. EFFECTIVE MOMENTUM APPROXIMATION 



Figure Q compares momentum distributions for the relativistic HS and NLSH models with the Woods- 
Saxon models used by Lapikas et al. 17] and by Dutta et al. The latter were based upon a NRDWIA 
analysis of the Saclay data for Tp = 99 MeV made by Bernheim et al. [2^. The differences between the 
HS and NLSH wave functions are relatively small, but the NIKHEF \p wave function has a much stronger 
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FIG. 1: (Color online) Momentum distribution for ^^C{e, e'p). 

peak and declines more rapidly with respect to Pm- Similarly, the NIKHEF Is wave function starts higher 
and falls faster than the relativistic wave functions. Interestingly, the Saclay wave functions show similar 
behavior but the effects, compared with the NIKHEF wave functions, are smaller for the Ip and larger for 
the Is wave function. The large spread among these momentum distributions will obviously be reflected in 
the corresponding DWIA calculations. 

To evaluate the effect of variations of p{pm) upon DWIA calculations using both relativistic and nonrel- 
ativistic wave functions, it is simplest to employ the EMA-noSV approach. Several previous papers have 
shown that the differences between EMA and full RDWIA calculations are relatively small for cross sections. 
Although that approximation does not provide the best description of the left-right asymmetry, Alt, it has 
the advantage that the nonrelativistic wave function can be used directly. 

Figure 121 compares EMA calculations using NLSH, HS, and NIKHEF wave functions with the JLab data 
for — 0.6, 1.2, and 1.8 (GeV/c)^. For each kinematics the top set are semi-inclusive data for Em < 80 
MeV, the second set shows a bin 15 < £",„ < 25 MeV dominated by the Ip shell, and the bottom set 
shows a bin 30 < Em < 50 MeV dominated by the Is shell. The curves were fit to data in the range 
\Pm\ < 200 MeV/c neglecting, at present, the mixing between the shells. We will later find that filling of the 
Ip minimum and broadening of the Is distribution can be described by mixin g. T hese calculations employ 
the EDADl optical potential; other choices of potentials from Cooper et al. [23| produce small variations 
in the normalization factors with little dependence upon the choice of overlap function and practically no 
discernible effect upon the curves. The results for the two relativistic wave functions, NLSH and HS, are 
very similar but the fitted pm distributions for the NIKHEF wave function tend to be too narrow, especially 
for the Ip shell. Consequently, the Ip normalization factors listed in Table ^ are systematically smaller 
for the NIKHEF wave function than for the NLSH or HS wave functions. These normalization factors are 
correlated well with the peak of the momentum distribution, with the NLSH results largest and NIKHEF 
results smallest. 

Similar comparisons are shown in Figs. I3I5I for several low experiments performed using quasiper- 
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FIG. 2: EMA-noSV calculations for ^^C{e,e'p) are compared with quasiperpendicular data for 0.6 < < 1.8 
(GeV/c)^ from JLab [isll . In each panel the top set of curves shows the inclusive 10 < Em < 80 MeV bin scaled by 
a factor of 20, the middle set shows the p-shell bin with 15 < Em < 25 MeV, and the bottom set shows the s-shell 
bin with 30 < Em < 50 MeV scaled by a factor of 0.1. Solid, dashed, and dash-dotted curves are based upon NLSH, 
HS, and NIKHEF wave functions, respectively. 

pendicular kinematics. Except when noted otherwise, a normahzation factor for each curve was obtained 
using a least-squares fit to the data for \pm\ < 200 MeV/c. Fig. O shows that the relativistic wave functions 
describe the data for the Saclay experiment with Tp — 87 MeV in quasiperpendicular kinematics better than 
the NIKHEF wave functions, which fall too rapidly with increasing The results shown in Fig. the 
Tokyo experiment with ^ o.29 (GeV/c)^ are somewhat ambiguous: although the NLSH and HS wave 
functions give better fits to for ^ 60 MeV/c, particularly for the ground-state, none of the calculations 
can describe the enhancement at lower pm seen for the ground state or the dip for the s-state. The low 
Pm bulge for the ground-state data seems somewhat implausible and it is likely that there were difficulties 
in cleanly separating the Ip and Is contributions to the low pm data. Therefore, the Ip normalization for 
this set was fit to 50 < Pm < 200 MeV/c. Similar calculations are also compared with the SLAG data for 

~ 1.1 (GeV/c)^ in Fig. We again find that the Pm distribution for the NIKHEF wave function is too 
narrow to fit the Ip data well. However, the Is data remain well above all of the calculations for large Pm- 
Lapikas et al. [l7| attributed this problem to \p contamination due to inadequate radiative unfolding. 

The normalization factors obtained for single-component EMA-noSV fits are listed in Table ^ improved 
two-component RDWIA fits will be presented in the next section. The Tokyo result is omitted because 
the Is admixture is appreciable. We find that these factors are correlated well with the peak value of the 
p{pm), such that the NIKHEF wave functions produce the smallest and the NLSH wave functions the largest 
normalization factors. Over a wide range of Q^, including both early data at low and recent data for 

~ 1 (GeV/c)^, the relativistic wave functions fit the p„i distributions relatively well while the NIKHEF 
wave functions are systematically narrower than the data. Thus, the data for quasiperpendicular kinematics 
consistently prefer the broader pm distributions and larger normalization factors obtained with Dirac-Hartree 
wave functions. 

The data for parallel kinematics, on the other hand, present special problems. We omit the Mainz data 
[S^ l because there appear to be significant normalization problems for that experiment The NIKHEF 
and Saclay data for parallel kinematics are compared with EMA-noSV calculations for NLSH, HS, and 
NIKHEF wave functions in Figs. [3and|Hl The results for both experiments are similar but none of these 
calculations reproduce the data for the Ip shell. Due to obvious the difficulties for pm < 0, the normalization 
factors listed in Table IIIII were fit to < Pm < 200 MeV/c. These factors are also correlated with the 
peak values of p{pm) but for the relativistic wave functions are somewhat lower than were obtained using 
quasiperpendicular kinematics. On the other hand, these fits underestimate the cross sections for pm < 
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FIG. 3: EMA-noSV calculations for ^^C(e,e'p) are compared with quasiperpendicular data for Tp — 87 MeV from 
Saclay |30|. Solid, dashed, and dash-dotted curves are based upon NLSH, HS, and NIKHEF wave functions, respec- 
tively. 
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FIG. 4; EMA-noSV calculations for ^^C(e,e'p) are compared with quasiperpendicular data for Tp = 99 MeV from 
Saclay |3Q). Solid, dashed, and dash-dotted curves are based upon NLSH, HS, and NIKHEF wave functions, respec- 
tively. 



and larger Q^. 

In order to improve standard NRDWIA calculations for parallel kinematics, van der Steenhoven et al. 
0,1111 resorted to adjusting the overlap function, the optical model, and the current operator simultaneously. 
The NIKHEF analysis used the nonrelativistic expansion of the current operator in powers oip/m proposed 
by McVoy and van Hove 34]; we verified that at second order this expansion, denoted NR2, gives results that 
are practically indistinguishable from our EMA-noSV calculations. Their analysis also used a nonrelativistic 
optical potential produced by Comfort and Karp [35|, denoted CK, and a Percy factor [s^ for the ejectile 
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FIG. 5: EMA-noSV calculations for ^■^C(c,e'p) are compared with quasiperpendicular data for Q'^ — 0.29 (GeV/c)'^ 
from Tokyo 




[MeV/c] 

FIG. 6: EMA-noSV calculations for ^^C(e,e'p) are compared with quasiperpendicular data for — 1.1 (GeV/c)^ 
from SLAG jSlJ. Solid, dashed, and dash-dotted curves are based upon NLSH, HS, and NIKHEF wave functions, 
respectively. 

with P = 0.85 fm. Figure 1^1 compares several variations of the NR2 calculations that are similar to those of 
van der Steenhoven et al. The dashed curve uses this optical potential and current operator with the upper 
component of the NLSH wave function and fails to reproduce the data for parallel kinematics with pm < 0; 
the results are similar to those shown in Fig. [7|in the EMA approach. The narrower NIKHEF momentum 
distribution improves the fit to the data but appears to be shifted to the left — a shift to the right was 
accomplished by modifying the optical potential. This modified potential, denoted MCO in Ref. jl^], was 
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TABLE II: Normalization factors for ^^C(e, e'p)(lp)~^ fit to quasiperpendicular data for the Ip bin using several 
overlap functions in the EMAnoSV approximation. Possible Is contributions were neglected. The first set uses the 
EDADl and the second the EDAIC optical potential. 



experiment 


Q2 

{GeY/cf 


Em 

MeV 


NLSH 


EDADl 
HS 


NIKHEF 


NLSH 


EDAIC 
HS 


NIKHEF 


Saclay 


0.16 


15-22 


0.637 


0.599 


0.483 


0.572 


0.542 


0.444 


Saclay 


0.18 


15-22 


0.757 


0.705 


0.554 


0.684 


0.640 


0.510 


JLab 


0.6 


15-25 


0.934 


0.878 


0.697 


0.909 


0.856 


0.682 


SLAC 


1.1 


15-25 


0.852 


0.791 


0.613 


0.887 


0.822 


0.633 


JLab 


1.2 


15-25 


0.904 


0.849 


0.671 


0.950 


0.890 


0.699 


JLab 


1.8 


15-25 


0.901 


0.854 


0.689 


0.886 


0.841 


0.679 



TABLE HI: Normalization factors for ^^C(e, e'p)(lp)^^ fit to data for parallel kinematics using several overlap func- 
tions in the EMAnoSV approximation. These results were extracted for the 15 < Em < 22 MeV bin using NLSH 
wave functions and the EDADl optical potential. 



experiment 


Tp 


NLSH 


HS 


NIKHEF 


NIKHEF 


70 


0.635 


0.616 


0.541 


Saclay 


99 


0.685 


0.636 


0.505 



intended to simulate channel coupling in the final state but our more detailed analysis of channel coupling 
[3^ gave much smaller effects for this reaction; hence, we regard this modification of the optical potential 
to be somewhat ad hoc. Finally, van der Steenhoven et al. inserted an enhancement factor for transverse 
components of the current operator and adjusted that factor to fit the peak for negative Pm- Applying this 
factor without further adjustment does provide a reasonably accurate fit to the reduced cross section for 
parallel kinematics. 

We have not attempted to tune this fit because we doubt that it is really possible to distinguish between 
variations of the overlap function, the current operator, and the final-state interactions using these data. 
Nor does it appear possible to fit the low data for both parallel and quasiperpendicular kinematics 
simultaneously using unique overlap functions and global optical potentials fit to proton elastic scattering. 
The relativistic wave functions provide goods fits to the Pm dependence of the data for quasiperpendicular 
kinematics even if the normalization factors for > 0.6 (GeV/c)^ are significantly larger than those for low 
Q^. On the other hand, the normalization factors for the two Saclay experiments with quasiperpendicular 
kinematics do not agree well with each other either. Therefore, the failure of EMA calculations with NLSH 
wave functions to reproduce the Pm dependence of data for parallel kinemtics with Tp < 100 MeV seems to 
be a more serious problem than the variation of normalization factors with Q^. Recognizing that off-shell 
ambiguities are more serious for parallel kinematics because x = Q'^/2muj varies over a wide range and 
expecting the reaction mechanism to be more reliable for large Q^, we believe that the reaction mechanism 
for parallel kinematics with Tp < 100 MeV is not sufficiently reliable to justify fitting the bound-state wave 
function and that it is likely that in such analyses variations of the overlap function would tend to compensate 
for deficiencies in the reaction model. 

If one wants to fit the momentum distribution for a discrete single-hole state to (e, e'p) data, it is probably 
best to use quasifree kinematics where (w, q) is constant and is as large as possible. The reliability of 
the reaction mechanism should improve as increases and the use of quasifree kinematics should minimize 
distortion of p{pm) due to variation of current operator off shell. It would then be of interest to test the 
reaction mechanism by comparison with data for parallel kinematics with large Tp. Ideally such tests should 
have sufficient resolution in Em to isolate discrete states without the complications of incoherent mixtures. 
Unfortunately, no high- resolution (e, e'p) data for parallel kinematics with Tp > 200 MeV and sufficient 
coverage in pm are available. 
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FIG. 7: EMA-noSV calculations for ^^C(e,e'p) are compared with data from from NIKHEF 18] using parallel 
kinematics for Tp = 70 MeV/c. 




FIG. 8: EMA-noSV calculations for ^^C(e, e'p) are compared with data from from Saclay using parallel kinematics 
for Tp = 99 MeV. Solid, dashed, and dash-dotted curves are based upon NLSH, HS, and NIKHEF wave functions, 
respectively. 
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FIG. 9: (Color online) NR2 calculations for ^^C(e, e'p) are compared with data from from NIKHEF [13 using parallel 
kinematics for Tp — 70 MeV/c. WS refers to the NIKHEF Woods-Saxon wave function and > shows the effect of 
enhancing transverse components of the current opertor. 
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FIG. 10: Contributions of IP3/2 and lsi/2 fitted to JLab data for ^^C{e,e'p). In each panel the top set of curves 
shows the inclusive 10 < Em < 80 MeV bin scaled by a factor of 20, the middle set shows the p-shell bin with 
15 < Em. < 25 MeV, and the bottom set shows the s-shell bin with 30 < Em < 50 MeV scaled by a factor of 0.1. The 
solid curves show the total fit, the dashed curves show the Ip contributions to the inclusive and p-shell bins, and the 
dash-dotted curves show the Is contributions to the inclusive and s-shell bins. In addition, scale-factor fits for the 
IPSM are shown as dotted curves for the inclusive bin but are practically indistinguishable from the solid curves. All 
calculations use NLSH wave functions and EDADl optical potentials in full RDWIA. 



IV. RDWIA ANALYSIS 



A. Normalization factors 



1. JLab 



Full RDWIA fits to the JLab data are shown in Fig. ^jusing NLSH wave functions, EDADl optical 
potentials, and the CC2 current operator in Coulomb gauge. For the inclusive bin, these fits are slightly 
better than the corresponding EMA-noSV fits shown in Fig. [3 because full RDWIA calculations reproduce 
the left-right asymmetry better than the EMA-noSV approximation 14]. These fits for the Ip and Is 
bins are significantly better than single-component fits because neither bin represents a pure single-particle 
configuration. The Is contribution to the Ip bin fills in the minimum and rounds the peaks. Similarly the 
Ip contribution to the Is bin broadens the momentum distribution — the reduced Is component provides 
a better fit for small Pm while the Ip component improves the fit to larger pm- However, for pm > 200 
MeV/c the data for both the s-shell and inclusive bins remain significantly above these two-component fits. 
Similar data for ^^0(e, e'p) [20ll38j | covering a broad range of missing energy show that for pm ^ 150 MeV/c 
there is much more strength at large Em than can be attributed to single-nucleon knockout from the s-shell. 
The continuum fraction increases with pm, which may be simulated in the two-component model by a Ip 
admixture. Assuming that the Is quasiholc momentum distribution is similar to the Dirac-Hartree wave 
function, it would appear that the continuum contribution to Em > 30 MeV has a broader momentum 
distribution than the 1^3/2 wave function used in these two-component fits. Although part of this may be 
contributed by Ip fragmentation, much of the fitted Ip strength is probably a surrogate for multinucleon 
continuum. Recent SCGF estimates suggest that approximately 68% of the correlated continuum is in partial 
waves with £ < 1 [S^ . No doubt higher partial waves begin to contribute to the continuum as Pm increases 
and it is likely that the £ = 1 continuum is also broader than the lp3 /2 quasihole momentum distribution. 

RDWIA fits to the JLab data using NLSH, HS, or NIKHEF wave functions are compared in Fig. 
The two Dirac-Hartree wave functions provide very similar fits with a slight preference for NLSH. Mixing 
between Ip and Is contributions improves fits to the Is bin using NIKHEF wave functions by flattening the 
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FIG. 11: RDWIA fits to the JLab data for JLab data for ^^C(e,e'p) are compared for NLSH (solid), HS (dasfied), 
and NIKHEF (dasii-dot) wave functions. In eacli panel the top set of curves shows the inclusive 10 < Em < 80 MeV 
bin scaled by a factor of 20, the middle set shows the p-shell bin with 15 < Em < 25 MeV, and the bottom set shows 
the s-shell bin with 30 < Em < 50 MeV scaled by a factor of 0.1. 

small pm and lifting the large Pm regions. However, two-component fits using NIKHEF wave functions were 
impossible for the Ip bin because the Is normalizations were negative; hence, the Ip NIKHEF fits shown here 
simplify scale RDWIA calculations for 1^3/2. Furthermore, the Ip NIKHEF fits are even worse and their 
normalization factors are even smaller in RDWIA than in EMA-noSV: with EDADl those normalization 
factors are 0.66, 0.63, and 0.66 for RDWIA compared with 0.70, 0.67, and 0.69 for EMA-noSV at the three 
settings. This reduction occurs because the RDWIA calculations for inherently nonrelativistic Woods- 
Saxon wave functions were made by replacing the Percy factor with the Darwin factor obtained using the 
spin-orbit potential; this procedure is explained in Ref. jl3j . However, because the spin-orbit potential 
used by van der Steenhoven et al. is stronger than that for NLSH or HS Dirac-Hartree wave functions, 
the smaller Darwin factor leads to a narrower momentum distribution for RDWIA than for EMA-noSV 
and, consequently, smaller normalization factors. The fits for the Ip bin overshoot the peaks and then fall 
too rapidly as Pm increases and these problems carry over, in slightly diluted form, to the inclusive bin. 
Therefore, we reject those parametrizations of the 1^3/2 and lsi/2 overlap functions. 

RDWIA normalization factors fit to the JLab reduced cross section data for ^^C{e,e'p) in three bins of 
missing energy are collected in Table HVl These results were obtained using EDADl optical potentials but 
the corresponding results for EDAD2, EDAD3, EDAIC, or EEI optical potentials differ by less than ±5% 
with systematics similar to Table im The differences between NLSH and HS results are small for 1^3/2 and 
modest for lsi/2- Most of the discussion will be based upon the NLSH results because they provide slightly 
better fits to the data. 

A rather large fraction, approximately 87%, of the IPSM Ip strength is found in the lowest Em bin. If we 
assume that fragmentation of the Ip-strength by collective modes is largely confined to excitation energies 
below 10 MeV, the 87% we find for the Ip-shell is consistent with SCGF estimates. The additional Ip 
strength fit to the reduced cross section for larger Em probably represents continuum contributions with 
i > that should not be added to the valence Ip strength. 

The Is contribution to the p-shell bin is more than predicted by Eqs. (|3I4|I . but could arise either from 

resolution, discrete ^'^ states, or a stronger low-Em tail. Approximately 65% of the Is IPSM strength 
appears in the 30 < Em < 50 MeV bin and corresponds to a net occupancy of 94% when scaled according to 
Eq. to account for the spreading width. However, this is probably an overestimate of the Is spectroscopic 
factor because one would expect approximately 10% of the Is strength to lie beyond 80 MeV missing energy; 
hence, if we scale Nis for the NLSH fit to the inclusive bin we would obtain more than 100%. Although 
the overestimation appears to be smaller for the HS fits, neither accounts for continuum contributions from 
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TABLE IV: Normalization factors for ^^C(e,e'p) fit to JLab data for three bins of missing energy and three values 
of Q2 in (GeV/c)^ The columns labeled ave is averaged with respect to Q^. The first set uses NLSH and the second 
HS wave functions. Both sets use EDADl optical potentials. 













NLSH 










HS 






state 


Em 


bin 


Q2 = 0.6 




1.2 


Q2 = 1.8 


ave 


= 0.6 




1.2 


Q2 = 1.8 


ave 


lP3/2 


15 - 


- 25 MeV 


0.885 


0.868 




0.847 


0.87 


0.852 


0.836 




0.830 


0.84 


30 - 


- 50 MeV 


0.060 


0.078 




0.077 


0.072 


0.063 


0.058 




0.062 


0.061 




10 - 


- 80 MeV 


1.021 


1.001 




0.977 


1.00 


1.008 


0.985 




0.973 


0.99 


lSl/2 


15 - 


- 25 MeV 


0.133 


0.102 




0.198 


0.14 


0.082 


0.050 




0.110 


0.08 


30 - 


- 50 MeV 


0.669 


0.616 




0.667 


0.65 


0.628 


0.575 




0.620 


0.61 




10 - 


- 80 MeV 


0.953 


0.887 




1.009 


0.95 


0.818 


0.757 




0.845 


0.81 


IPSM 


10 - 


- 80 MeV 


0.981 


0.991 




0.984 


0.985 


0.949 


0.914 




0.940 


0.934 



rescattering or two-body currents that could artificially enhance Ni^. The difficulties in reproducing the cross 
sections for pm > 200 MeV/c in the E„i > 30 MeV region also suggest additional continuum contributions 
with broader momentum distribution that are not described by these two-component fits. Therefore, it is 
reasonable to assume that the depletion of the Is shell by short-range correlations is similar to that for the 
Ip shell and to attribute the additional strength in the inclusive cross section to continuum processes that 
are not related directly to the single-hole spectral function. 

The final fine of Table Hvl lists IPSM scale factors fit to the Em < 80 MeV data according to Eq. lO; 
these results are in reasonable agreement with the weighted average of the Ip and Is normalization factors. 
The corresponding fits arc shown in Fig. 1101 as dotted curves that are practically indistinguishable from 
the two-component fits to the inclusive data. Therefore, we find approximately 98.5% of the IPSM strength 
in Era < 80 MeV, which is considerably more than the 85% one expects after accounting for short-range 
correlations. It is likely that much of the extra strength observed for E^ > 30 MeV arises from processes 
other than single-nucleon knockout from the Is-shell, such that the total single-hole strength for Em < 80 
MeV is overestimated by this analysis. 

Under these conditions, we consider A^ip(15, 25) — 0.87 to be the most reliable estimate of the depletion of 
IPSM orbitals in ^^C by short -range correlations. This estimate is consistent with the direct measurements 
of correlated continuum made by Rohe et al. ;4| using parallel kinematics and the analysis of those data by 
Miither and Sick J39J based upon SCGF. The extra yield observed by Dutta et al. |15iJ in quasiperpendicular 
kinematics should then be attributed to two-body currents and rescattering processes in which some of the 
flux absorbed by the optical potential is distributed to final states that remain within the experimental 
acceptance. Several studies suggest that rescattering contributions to the continuum are minimized for par- 
allel kinematics 0, l40l |4]| . These background contributions affect the interpretation of nuclear transparency 
measurements, which will be considered in Sec. IIVBI 

For the sake of completeness, we also show in Fig. ^|thc left-right asymmetry in reduced cross section 

ffrodC^ = 0) - (T^cdid) = Tt) 
O-rod(0 = 0) + 0-rod((?> = TTj 

where the azimuthal angle = corresponds to an ejectile momentum in the electron scattering plane 
between the beam direction and the momentum transfer. The sensitivity of this quantity to spinor distortion 
was studied in some detailed in Ref. 0| for pure single-particle configurations. Figure [T^ shows that the 
lP3/2 admixture for the s-shell bin provides slightly better fits to the data for pm ^ 300 MeV/c. That effect 
is small for the p-shell bin because the Ip admixture is small and is also small for the inclusive bin because 
the asymmetries for the two components are similar, depending more upon the Dirac potentials than upon 
the details of the bound-state wave functions. Two-component fits based upon EMA-noSV calculations 
demonstrate that ~ without spinor distortion because the characteristic left-right asymmetry for 
electron scattering by a moving free proton is removed by using the reduced cross section. The remaining 
asymmetry due to dynamical enhancement of lower components of Dirac spinors is described very well by 
RDWIA calculations. As argued in Ref. p^, the flattening of the pm distribution for a^T in the s-shell 
and inclusive bins for increasing is probably caused by continuum contributions that need not share the 
characteristic asymmetry of single-nucleon knockout. 
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FIG. 12: RDWIA calculations for the left-right asymmetry in reduced cross section for quasiperpendicular kinematics 
with 0.6 < < 1.8 {GeV/cf. NLSH wave functions, EDADl optical potentials, and the CC2 current operator were 
used. The dashed curves show single-component calculations while results for two-component Ip + Is fits are shown 
as solid curves. The dash-dotted curves show two-component fits in EMA-noSV. 



2. Previous experiments 

The RDWIA analysis of the JLab data presented in the preceding section provides a satisfying degree of 
internal consistency over a broad range of Q^. Unfortunately, similar analyses for older data sets present 
a variety of consistency problems that cannot be resolved by simply by choosing the best overlap functions 
or optical potentials. Therefore, in this section we consider those data in chronological order using only a 
single set of options, namely NLSH and EDADl, in the RDWIA framework. 

The fits shown in Fig. ^|appear to describe the Tokyo data (2^ very well, but the Is contribution to the p- 
shell bin is implausibly large. With NLSH wave functions and EDADl optical potential, we find {Nip, Nis) = 
(0.56, 0.42) for 6 < E,n < 30 MeV and {Nip, Ni^) = (0.08, 0.69) for 21 < E„, < 66 MeV. Apparently there was 
an experimental problem in defining the p-shell or there is an unidentified problem in the reaction mechanism 
that enhances the cross section for small Pm- Because this is the only quasiperpendicular data set with severe 
problems in the p-shell bin for small Pm-, we believe that energy resolution is a more likely explanation than 
reaction mechanism. The single-component fit to Pm > 50 MeV shown in Fig. [S] corresponds to Nip = 0.77, 
which is more consistent with the other quasiperpendicular data for this bin. The Ip contribution to the Is 
bin, on the other hand, is fairly small. Scaling the Is fit for 21 < E„i < 66 MeV according to Eq. gives 
a net occupancy of 80%, which is reasonable but smaller than the corresponding JLab result. Perhaps the 
continuum contamination was smaller for the Tokyo experiment. 

RDWIA provides no qualitative improvement over EMA-noSV fits to NIKHEF or Saclay data for parallel 
kinematics. As discussed before, we believe the severe variation of Bjorken x in parallel kinematics distorts 
momentum distributions fit to (e, e'p) data for low using familiar off-shell extrapolations of the single- 
nucleon current operator. Therefore, we do not consider those data further. Conversely, single-component 
RDWIA fits to the Saclay data for quasiperpendicular kinematics are similar to those shown in Figs. I3I4I 
There is no evidence for Is contamination of the 15 < Em < 25 MeV bin and the 30 < Em < 50 MeV bin 
is too narrow and the range of Pm is too small to find significant Ip strength there. Furthermore, mixed fits 
sometimes produce negative contamination factors with large uncertainties; hence, only single-component 
fits are considered for the Saclay quasiperpendicular data. We find Ip normalization factors of 0.63 for Ref. 
[30l | or 0.75 for Ref. |2^ that are somewhat smaller than the JLab results but are not in especially good 
agreement with each other either. Similarly, Nis = 0.46 for 30-50 MeV corresponds to an occupancy of 
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FIG. 13: RDWIA fit to Tokyo data |29| using NLSH wave functions. Solid curves show the fits while dashed and 
dash-dotted curves show Ip and Is contributions. 



67%, which is also somewhat less than the JLab result. Given that we believe the reaction mechanism is 
more reliable at higher and that there are unresolved problems in reproducing data in parallel kinematics 
at the same Q^, we are not especially alarmed by a 10-20% discrepancy between < 0.2 and > 0.6 
(GeV/c)2. 

Finally, the statistical quality of the Ip SLAG data for = 1.1 (GeV/c)^ do not support an extraction 
of a Is admixture either. Nor did we attempt to fit the Ip contribution to the 30 < < 80 MeV bin 
because the enhancement for large positive Pm is probably an artifact of inadequate radiative unfolding. 
Single-component fits using RDWIA calculations are very similar to the EMA fits shown in Fig. Eland the 
figures are omitted. The SLAG Ip normalization factor for the p-shell bin is 0.86 for NLSH and EDADl, 
which is perfectly consistent with the JLab results. However, Nis — 1.2 for 30-80 MeV is much larger 
than the JLab result and scaled up to unrestricted Em would represent 140% of IPSM, suggesting a rather 
substantial underlying continuum. 



3. dependence of spectroscopic factors 

Lapikas et al. proposed that the variation seen in their analysis of the spectroscopic factors for ^^G 
could be explained in terms of the resolution with which the electromagnetic knockout reaction probes the 
structure of a quasiparticle, such that the true spectroscopic factor would be observed in the limit — > oo. 
Although this model was expressed in terms of a scale-dependent renormalization of the spectroscopic factor, 
it would be more appropriate to attribute such an effect, if present, to the single-nucleon current operator 
instead of to the spectroscopic factor. The overlap function for single-nucleon knockout takes the form 

(*i-'-^^l*^^^)=E^"/30,(r) (7) 

where is a many-body state of either the A or A — 1 nucleus, a is a state label with being the target 
ground state, is a parentage coefficient, and t/)^ is a local overlap function that is expected to resemble 
an orbital for the mean field. If the expansion basis is orthnormal, the spectroscopic factors become 



S(3 = ^ |Cq/3| 
a 



(8) 
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The parentage expansion depends upon the nuclear structure and is independent of the probe; hence, the 
spectroscopic factors should not depend upon . However, the assumption that the electromagnetic current 
operator is adequately approximated by the free nucleon current with the energy transfer placed on shell and 
with free form factors is questionable. Such approximations violate current conservation, Lorentz covariance, 
and unitarity. Furthermore, the current in a many-body system becomes nonlocal when expressed in terms 
of single-nucleon degrees of freedom; such nonlocality could contribute to an apparent scale dependence. 
Nevertheless, one must resort to simplified models of the current operator because there is no practical 
method for performing rigorous calculations for nucleon knockout from nuclei with A > 2. Therefore, 
dependent modifications of the off-shell single-nucleon current operator should not be surprising. 

The recent JLab data for ^^C show no evidence for systematic variations of the current operator within 
the range 0.6 < Q2 < 1.8 (GeV/c)2. Although consistency among previous experiments is not entirely 
satisfactory, our RDWIA analysis shows much less variation with than suggested by Lapikas et al. 
Because the continuum contamination for Em > 30 MeV has not been modeled accurately, we consider Nip 
for the p-shell bin to be the most reliable gauge of possible dependencies. Comparing the average of 
the two Saclay results for quasiperpendicular kinematics with the average of the JLab results for NLSH and 
EDADl, we see an increase of 0.87/0.69 ~ 1.26 compared with a factor of about 1.6 for Sip + Sis from Ref. 
pjj. There may still be a small effect, which we would attribute to the current instead of spectroscopic 
factors, but until the failure of RDWIA to reproduce low data for parallel kinemtics is understood we 
are reluctant to interpret this as a simple scale factor. 



B. Transparency 

The experimental definition of nuclear transparency 

/ /y d^PmdEmdile K Paip„i)Sa{E„i)aep 

compares the measured semi-inclusive differential cross section for proton knockout with a PWIA calculation 
based upon factorization of a model spectral function and the free e-p cross section. The phase-space volume 
is indicated by V and we leave possible correction factors for nonuniform acceptance implicit, assuming that 
they are applied in both numerator and denominator consistently. Here K is a, kinematic factor, Pa{Pm) 
is the momentum distribution for orbital a, Sa{Em) is the energy distribution for single-nucleon removal 
from orbital a, and / is a correction factor describing the depletion of the single-hole spectral function by 
correlations that shift strength to large Em- Dutta et al. assumed that / is independent of a and obtained 
an estimate of 0.9 for ^^C from 31] that is consistent with the Ip normalization factor we extracted in the 
preceding section. 

Similarly, a theoretical definition of nuclear transparency appropriate to quasiperpendicular kinematics 
takes the form 

j-^ _ J2a I dpmPm CTRpwi A (Pm , -Eg ) ^^^-^ 
J2al dprnPm 0-RPWIA (Pm , ) 

where the numerator integrates the RDWIA and the denominator the RPWIA cross section for the same 
model of the spectral function. For simplicity we neglect the spread of electron kinematics and assume 
quasifree kinematics for the ejectile, neglecting the small variation of transparency with nucleon kinetic 
energy within the acceptance. Similarly, we also neglect variations over the spreading widths of nuclear 
orbitals and assume that the bound-nucleon kinematics can be approximated by the IPSM. A more rigorous 
calculation would use Monte Carlo integration over a realistic model of the experimental acceptance, but the 
dependence of nuclear transparency upon energy is very mild and the differences between transparencies for 
parallel and quasiperpendicular kinematics are small |42| |. 

Note that one should not include the correlation factor / in the theoretical definition of transparency be- 
cause any strength that is removed from the numerator by correlations is also removed from the denominator 
by using the same spectral function for both. This factor is appropriate for Eq. (PJ because the numer- 
ator integrates experimental yield while the denominator integrates a theoretical calculation based upon a 
model spectral function. When the denominator uses the IPSM one requires a theoretical estimate for the 
fraction of the the spectral strength that is shifted out of the experimental range of Em- Unfortunately, a 
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FIG. 14: RDWIA calculations using several optical potentials are compared to transparency data for ^^C(e, e'p). The 
solid curve has been multiplied by 1.12 to provide a better fit to the data from MIT (open circle), SLAG (open square) 
and JLab (filled circles. The SLAG datum includes systematic uncertainties, but other error bars are statistical only. 



recent comparison of the data with RDWIA and Glauber calculations includes the correlation factor in both 
theoretical and experimental definitions of transparency |4^ . which we believe is incorrect. 

RDWIA calculations using NLSH wave functions and several optical potentials from Dirac phenomenology 
[2^ are compared in Fig. El with transparency data for -^^C(e, e'p) from MIT El. S LAG [Sllli^, and JLab 
|l5l l46l | . The present results are very similar to those we previously obtained jl^l using the EMA-noSV 
approximation and a definition of transparency that replaces the differential cross sections in Eqs. (I9I10|) 
with the corresponding reduced cross sections. Neither change is more than a couple percent. Our results 
are also similar to the RDWIA results of Lava et al. 43] after removal of the correlation factor. Also note 
that Lava et al. used a transparency definition based upon reduced cross section. However, Meucci [4^ 
obtained somewhat larger transparency factors using an EMA-SV calculation based upon differential cross 
section. The latter are the only RDWIA calculations, of which we are aware, that slightly overestimate the 
transparency for ^^G{e, e'p), but the origin of their enhancement is not known. 

The variation among >l-independent global optical potentials is small, but the yl-dependent potential fit 
to data for ^'^G{p,p) provides the best description of the energy dependence for Tp < 400 MeV. However, all 
of these calculations systematically underestimate the experimental transparency. Multiplying the EDAIC 
results by 1.12 provides a reasonable description of the data over this range of ejectile energy. Interestingly, 
this factor is very close to the inverse of the correlation correction applied by Abbott et al. If we 

assume that their correction factor is a reasonable estimate of the fraction of the single-nucleon yield shifted 
by correlations to large E^, this comparison suggests that the experimental yield is approximately 12% larger 
than single-nucleon knockout. This is consistent with the fact that the RDWIA normalization factors sum to 
the full IPSM strength without depletion. Therefore, either the RDWIA is approximately 12% too absorptive 
or there is a significant continuum arising from multinucleon knockout that contributes approximately 12% of 
the strength for E„i < 80 MeV. The rather strong ^*^0(e, e'p) continuum seen in Fig. 21 of Ref. suggests 
that the underlying continuum is probably responsible for the large normalization factors and transparency 
values for ^^C{e,e'p). 

The consistency between the fitted Ip normalization factors for > 0.6 (GeV/c)^ and the predicted 
IPSM depletion factor supports the accuracy of (e, e'p) attenuation factors calculated using RDWIA with 
global optical potentials. The observation that the yield for 30 < Em < 80 MeV is considerably larger 
than can be attributed to Is knockout demonstrates that there is significant background in the numerator 
that is not described by the spectral function used in the denominator of Eq. The contribution 
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to this background from multinuclcon currents artificially increases the measured transparency insofar as 
that quantity is intended to represent the loss of flux from single-nucleon knockout channel. However, 
the contribution to this background from rescattering processes that simply redistribute that flux within 
the experimental acceptance, instead of removing it entirely, represents a limitation of the distorted wave 
approximation which does not account for where flux "absorbed" by the optical potential finally appears 
|42|. The background fraction probably increases with A, but a more detailed model of the continuum is 
needed to distinguish between the processes that increase the experimental transparency with respect to 
RDWIA. 



V. SUMMARY AND CONCLUSIONS 

We have used RDWIA calculations based upon Dirac-Hartree bound-state wave functions to analyze data 
for i2c(e,e'p) for < 80 MeV. Good fits are obtained for the recent JLab data (Tsl for 0.6 < < 1.8 
(GeV/c)^, with a slight preference for NLSH over HS wave functions and practically no variation with any 
of the optical potentials provided by Cooper et al. . We find that the p-shell bin contains approximately 
87% of the IPSM strength independent of over that range. Although the s-shell bin appears to carry 
almost 100% of IPSM, independent of Q^, there is evidence in the momentum distribution for a significant 
continuum that artificially increases the Is normalization factor for two-component fits mixing f^^^d 
lsi/2 contributions. Therefore, we consider the Ip contribution to the 15 < Em < 25 MeV bin to be the 
most reliable estimate of the depletion of IPSM orbitals by short-range correlations; the occupancy of 87% 
is consistent with estimates based upon the self-consistent Green's function method and with recent direct 
measurements of the correlated continuum using parallel kinematics 0, ^9. 

We have also analyzed the same low data sets used by Lapikas et al. jl^l to study the dependence of 
normalization factors for ^^C(e, e'p). We find that the RDWIA calculations for Dirac-Hartree wave functions 
reproduce the low data for quasiperpendicular kinematics well, but with somewhat smaller normalization 
factors. For example, the Ip normalization factors for two Saclay experiments were 0.63 and 0.75 using 
NLSH wave functions and EDADl optical potentials, compared with 0.87 for the JLab data. However, the 
same model fails to reproduce low data for parallel kinematics and remediation of this problem will 
require more than a multiplicative factor. Similar problems in previous NRDWIA analyses prompted van 
der Steenhoven et al. |l8 | to adjust the Woods-Saxon wave function, the optical potential, and an empirical 
enhancement of the transverse components of the current simultaneously. Although a better fit was achieved, 
this wave function does not reproduce the data for quasiperpendicular kinematics for either early low 
experiments or the more recent high experiments. The momentum distributions for the fitted Woods- 
Saxon wave functions are too narrow, which artificially reduces the normalization factors fit to data. Thus, 
Lapikas et al. ^3 obtained a Ip normalization factor of only 0.56 for < 0.2 (GeV/c)^ and suggested that 
spectroscopic factors might depend strongly upon the resolution of the probe. We argue that spectroscopic 
factors are properties of nuclear structure that are independent of probe, but that the effective current 
operator for single-nucleon knockout may include additional dependencies beyond those in familiar off- 
shell extrapolations of the free nucleon current operator. On the other hand, the additional dependence 
is probably smaller than their estimate. Further investigation is needed to explain the discrepancy between 
parallel and quasiperpendicular kinematics at low — we believe that it is unwise to fit the momentum 
distribution to data with a strong variation of x = (5^/2maj in the absence of a more fundamental theory of 
the effective current operator. 

Finally, we found that RDWIA calculations of nuclear transparency need to be multiplied by approxi- 
mately 1.12 to reproduce data for ^^C and attribute much of this enhanced transparency to contributions 
to the continuum for 30 < E,n < 80 MeV that are not directly related to single-nucleon knockout, such 
as multinuclcon electromagnetic currents. However, a more detailed model of the continuum is need to 
distinguish between contributions due to multinuclcon currents and rescattering processes. 
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